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ABSTRACT. The behavior at criticality of spatial SIR (susceptible/infected/recovered) epidemic mod- 
els in dimensions two and three is investigated. In these models, finite populations of size N are situ- 
ated at the vertices of the integer lattice, and infectious contacts are limited to individuals at the same 
or at neighboring sites. Susceptible individuals, once infected, remain contagious for one unit of time 
and then recover, after which they are immune to further infection. It is shown that the measure- 
valued processes associated with these epidemics, suitably scaled, converge, in the large- AT limit, 
either to a standard Dawson- Watanabe process (super-Brownian motion) or to a Dawson- Watanabe 
process with location-dependent killing, depending on the size of the the initially infected set. A key 
element of the argument is a proof of Adler's 1993 conjecture that the local time processes associated 
with branching random walks converge to the local time density process associated with the limiting 
super-Brownian motion. 

1. Introduction 

1.1. Spatial SIR epidemics. Simple spatial models of epidemics are known to exhibit critical 
thresholds in one dimension: Roughly, when the density of the initially infected set exceeds a certain 



level, the epidemic evolves in a m arkedly different fashion than its b ranc hing envelope. S ee Lalley 



2007 ) for a precise statement, and lAldousI ( 1997h jMartin-Loll ( 1998h . and lDolgoarshinnykh and Lalley 



20061) for analogous results in the simpler setting of mean-field models. The main purpose of this 



article is to show that spatial SIR epidemics (SIR stands for susceptible/infected/recovered) in di- 
mensions two and three also exhibit critical thresholds. 

The epidemic models studied here take place in populations of size N located at the sites of the 
integer lattice Z d in d dimensions. Each of the N individuals at a site x 6 Z d may at any time be 
either susceptible, infected, or recovered. Infected individuals remain infected for one unit of time, 
and then recover, after which they are immune to further infection. The rules of infections are as 
follows: at each time t = 0, 1, 2, . . . , for each pair (i x , s y ) of an infected individual located at x and 
a susceptible individual at y, i x infects s y with probability pn{x, y)- We shall only consider the case 
where the transmission probabilities pw{x,y) are spatially homogeneous, nearest-neighbor, and 
symmetric, and scale with the village size N in such a way that the expected number of infections 
by a contagious individual in an otherwise healthy population is 1 (so that the epidemic is critical), 
that is, 

Assumption 1. pn(x; y) = l/[(2d + l)N] if \y — x\ < 1; and = otherwise. 

Our main result, Theorem [2] below, asserts that under suitable hypotheses on the initial configu- 
rations of infected individuals, the critical spatial SIR—d epidemic can be rescaled so as to converge 
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to a Dawson-Watanabe measure-valued diffusion in both d = 2 and d = 3. Depending on the size 
of the initially infected set, the limiting Dawson-Watanabe process h as eithe r a pos itive killing rate 



or no killing at all. The analogous result for d = 1 was proved in iLalleyl (120071) . using the fact 
that one-dimensional super-Brownian motion (the Dawson-Watanabe process with no killing) has 
sample paths in the space of absolutely continuous measures. In higher dimensions this is no longer 
true, so a different strategy is needed. 

1.2. Branching envelope of a spatial epidemic. The spatial SIR epidemic in d dimensions is natu- 
rally coupled with a nearest neighbor branching random walk on the integer lattice 7L d \ this branch- 
ing random walk is often referred to as the branching envelope of the epidemic. Particles of this 
branching random walk represent infection attempts in the coupled epidemic, some of which may 
fail to be realized in the epidemic because the targets of the attempts are either recovered or are 
targets of other simultaneous infection attempts. The branching envelope evolves as follows: Any 
particle located at site x at time t lives for one unit of time and then reproduces, placing random 
numbers £ y of offspring at the sites y such that \y — x\ < 1. The random variables £ y are i.i.d., 
with Binomial(iV, l/[(2d + l)iV]) distributions. Denote this reproduction rule by TZn, and denote 
by IZoo the corresponding offspring law in which the Binomial distribution is replaced by the Pois- 
son distribution with mean 1/ (2d + 1). Since offspring are placed independently at each of the 
2d+ 1 nearest neighbors, the expected total number of offspring of a particle is 1, i.e. the branching 
random walk is critical. Moreover, under either reproduction rule 1Zn or TZoo, particle motion is 
governed by the law of the simple nearest neighbor random walkQ on 7L d (with holding probability 
l/(2d + 1)): In particular, given that a particle at site x has k offspring, each of these offspring 
independently chooses a neighboring site y according to the law 

(1) Pi(x,y) = l/{2d + 1) for |y-x|<l. 

Note that the covariance matrix of the increment has determinant a 2d , where a 2 is the variance 
parameter of the jump distribution, defined by 

o ( 2 

(2) a 2 :-- 



2d+ 1 

The spatial SIR-d epidemic can be constructed together with its branching envelope on a common 
probability space in such a way that the branching envelope dominates the epidemic, that is, for each 
time n and each site x the number of infected individuals at site x at time n is no larger t han the 



numbe r of particles in the branching envelope. The construction, in brief, is as follows (see lLalley 



(120071) ): Particles of the branching random walk will be colored either red or blue according to 
whether or not they represent infections that actually take place, with red particles representing 
actual infections. Initially, all particles are red. At each time t = 0,1, 2,..., particles produce 
offspring at the same or neighboring sites according to the law TZn described above. Offspring 
of blue particles are always blue, but offspring of red particles may be either red or blue, with the 
choices made according to the following procedure: All offspring of red particles at a location y 
choose numbers j G [N] := {1,2, .. . ,N} at random, independently of all other particles. If a 
particle chooses a number j that was previously chosen by a particle of an earlier generation at the 
same site y, then it is assigned color blue. If k > 1 offspring of red particles choose the same number 



'Throughout the paper, the term simple random walk will mean simple random walk with holding probability l/(2d + 

1). 
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j at the same time, and if j was not chosen in an earlier generation, then 1 of the particles is assigned 
color red, while the remaining k — 1 are assigned color blue. Under this rule, the subpopulation of 
red particles evolves as an SIR—d epidemic. 

It is apparent that when the numbers of infected and recovered individuals at a site and its nearest 
neighbors are small compared to N, then blue particles will be produced only infrequently, and 
so the epidemic process will closely track its branching envelope. Only when the sizes of the 
recovered and infected sets reach certain critical thresholds will blue particles start to be produced 
in large numbers, at which point the epidemic will begin to diverge significantly from the branching 
envelope. Our main result, Theorem [2] below, implies that the critical threshold for the number of 
initially infected individuals is on the order N 1 ^ 3 ^ d ^ 2 \ 

The SIR-d epidemic is related to its branching envelope in a second — and for our purposes more 
important - way. The law of the epidemic, as a probability measure on the space of possible popu- 
lation trajectories, is absolutely continuous relative to the law of its branching envelope. The like- 
lihood ratio can be expressed as a product over time and space, with each site/neighbor/generation 
contributing a factor (see § !3.3l below). Each such factor involves the total occupation time (x) 
of the site, that is, the sum of the number of particles at site x over all times prior to n. Thus, the 
asymptotic behavior of the occupation time statistics for branching random walks will play a central 
role in the analysis of the large- N behavior of the SIR-d epidemic. 



1.3. Watanabe's Theorem. A fundamental theorem of Watanabe (1968) asserts that, under suit 



able rescaling (the Feller scaling) the measure-valued processes naturally associated with critical 
branching random walks converge to a limit, the standard Dawson-Watanabe process, also known 
as super-Brownian motion. 

Definition 1. The Feller-Watanabe scaling operator J-f. scales mass by 1/k and space by 1/Vk, 
that is, for any finite Borel measure fi(dx) on M. d and any test function tp £ C^°(M. d ), 

(3) F h n) = AT 1 / ^(Vkx)n(dx). 



Watanabe's Theorem . Fix N, and for each k = 1, 2, . . . , let X k be a branching random walk 
with offspring distribution IZn and initial particle configuration Xq. (In particular, X k (x) denotes 
the number of particles at site x G Z rf in generation [t], and X k is the corresponding counting 
measure.) If the initial mass distributions converge, after rescaling, as k — > oo, that is, if 

(4) F k X k ^p = X 

for some finite Borel measure p, on W 1 , then the rescaled measure-valued processes {J~kX^)}.t 
converge in law as k — > oo: 

{FkX k )kt x t , 

where =>■ represents the weak convergence relative to the Skorokhod topology on 

D([0, oo); Mp{E. d )). The limit is the standard Dawson-Watanabe process X t (super-Brownian 

motion) with variance parameter a 1 (equivalently, standard super-Brownian motion run at speed 

a). 



See lEtheridge d2000h for an in-depth study of the Dawson-Watanabe process and a detailed proof 



of Watanabe's Theorem. Because the process X t has continuous sample paths in the space of finite 
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Borel measures, it follows routinely from Watanabe's theorem that the occupation measures for 
branching random walks converge to those of super-Brownian motion: 



Lemma 1. The following joint convergence holds: 



(5) 



X, 



X s ds 



{F k X K ) kt ,( ( (F k X k ) ks ds 

\Jo // \ JO 

where represents weak convergence relative to the Skorokhod topology on D([0, oo); Mp(M. d )) 2 . 

Proof. The Dawson-Watanabe process X t has continuous sample paths in D([0, oo); Mp(W i )), see, 
e.g., Proposition 2.15 in lEtheridgel(l200nh . The functional (X t ) i-> (/* X s ds) is continuous relative 
to the Skorokhod topology on the subspace of continuous measure-valued processes, so the result 
follows from Watanabe's theorem and the continuous mapping principle. □ 

1.4. Local times of critical branching random walks. In dimension d = 1 the super-Brownian 
motion has sample paths in the space of absolutely continuous measures, that i s, for each t > 
the random mea s ure X t is absolutely continuous relative to Lebesgue measure ( Konno and Shigal 
( 1988 ). Reimers (198^)). More over, the Rado n-Nikodym derivative X(t, x) is jointly continuous in 
t, x (for t > 0). It is shown in lLalleyl (|2007f) that if a sequence of branching random walks satisfy 
the assumptions in Watanabe's Theorem, then the density processes associated with those branching 
random walks, under some smoothness assumptions on the initial configurations and after suitable 

scaling, converge to the density process of the limiting super-Brow nian motion. 

In dimensions d > 2 the measure X t is almost surely singular ( Dawson and Hochberd ( 1979h ). 
Therefore, one cannot expect the convergence of density processes as in lLalleyl (I2007h . We shall 
prove, however, that the occupation measures of critical branching random walks have discrete 
densities that converge weakly — see Theorem \T\ below. The limit process is the local time process 
associated with the occupation measure 



Li 



I", ds 



of the super-Brownian motion. In dimensions d 
absolutely continuo us, despite the fact that X t is singular 



2, 3, the rand om measure L f i s, for each t > 0, 
Susitanil d 19891) . Ilscod d 19861) and 



see 



Fleischmannl dl988h . Moreover, under suitable hypotheses on the initial condition Xq, the density 
process L t (x) is jointly continuous for t > and x G U. d : This is the content of Sugitani's theorem. 
For the reader's convenience, we state Sugitani's Theorem precisely here. For t > and x £ M. d , 
set 

ft e -M 2 /2i 
q t {x) = / (f) s (x)ds, where cp t (x) 



(27Tt) d / 2 



is the usual heat kernel. 



Sugitani's Theorem . Assume that d = 2 or 3, and that the initial configuration /i := Xq of the 
super-Brownian motion Xt is such that the convolution 

(6) (qt * fJ-)(x) is jointly continuous in t > and x € M. d . 

Then for each t > the occupation measure Lt is absolutely continuous, and there is a jointly 
continuous version Lt(x) of the density process. 
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We call (L t (x)) t>Q x&Rd the local time density process associated with the super-Brownian mo- 
tion. In view of Watanabe's and Sugitani's theorems, it is natural to conjecture (see Remark Q] 
below) that the local time density processes of branching random walks, suitably scaled, converge 
to the local time density process of the super-Brownian motion. Theorem Q] below asserts that this 
conjecture is true. Let X k be a sequence of branching random walks on Z d . Write 

(7) X*(x) : = # particles at x at time i, and 

**(*): = £jtf(x). 

(We use the notation instead of L\ because in the corresponding spatial epidemic model, the 
quantity R k (x) represents the number of recovered individuals at site x and time n.) Denote by 

the transition probability kernel of the simple random walk on Z d , that is, P n = P * p n_1 is the nth 
convolution power of the one-step transition probability kernel given by ([]]). Let G n {x,y) be the 
associated Green's function: 

G n {x) :=Y,Pi{x). 

For any finite measure u on Z rf with finite support, set 

(fiG n )(x) := (n* G n )(x) = ^2 fi(y)G n (x - y), 

y 

and denote by ^Gt{y) the continuous extension to [0, oo) x M, d by linear interpolation. 

Theorem 1. Assume that d = 2 or d = 3. For each k = 1, 2, . . . , let X k be a branching random 
walk whose offspring distribution is Poisson with mean 1. Assume that the initial configurations 
fi k := Xq satisfy hypothesis (HJ) of Watanabe's theorem, where the limit measure u has compact 
support and satisfies the hypothesis © of Sugitani's theorem. Assume further that 

ra\ fi k G kt (Vkx) 2 

(9) — — => iw<x 2 t * M/v \( x )' 

where indicates weak convergence in the topology of D([0, oo), Cfo(M rf )). Then as k — > oo, 

where Lt(x) is the local time density process associated with the super-Brownian motion with vari- 
ance parameter a 2 started in the initial configuration Xq = \x. 

Theorem Q] will be proved in §2. 



Remar k 1. The analogous result for critical branching Brownian motions was conjectured by lAdler 



(119931) . who proved the marginal convergence for any fixed t and x. 



Remark 2. The assumption that the offspring distribution is Poisson with mean 1 can be relaxed. 
All that is really needed is that the offspring distribution has an exponentially decaying tail. See 
Remark |7J 
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Remark 3. The hypothesis (0]) does not by itself imply (©, even if the limit measure fi satisfies the 
hypothesis © of Sugitani's theorem. Sufficient conditions for Q are given in Proposition [Qbelow. 
In particular, in dimension 2, if (O holds and the maximal number of particles on a single site is 
bounded in k, then © is satisfied. 

Remark 4. Let X t be super-Brownian motion in dimension d = 2. For each t > the random 
measure X t is singular, so by Fubini's theorem, for almost every point x € M? the set of times 
t > such that x is a point of density of X t has Lebesgue measure 0. Under hypothesis © we can 
make an analogous quantitative statement for branching random walk: For any fixed x G 1? and all 

t > 0, 

[kt] 

(11) E^I x ^ x)>0 = O(k/logk). 

m=l 



Proof. By Proposition 35 in lLalley and Zhengl (|2007r) . there exists 5 > such that for all k and m 
sufficiently large, 

~X k m (x)\X k m (x) > 



E 



> dlogm. 



But hypothesis Q implies that 



ER k kt {x) = { l x k G kt ){x)=0{k) 1 



and 



ER k t {x) = ]T EX k m {x) =Y, E [x k m (x)\X k m {x) > 0] • P [x k m {x) > 



m<kt 



m<kt 



□ 



1.5. Scaling limit of spatial SIR epidemic. Before stating our result, we first recall the definition 
of Dawson- Watanabe processes with variable-rate killing. The Dawson-Watanabe process X t with 
killing rate 9 = 9(x, t, m) (assumed to be progressively measurable and jointly continuous in (t, x)) 
and v ariance parameter a 1 can be characterized by a martingale problem (IDawson and Perkins 
dl999h . §6.2): For any test function ^eC c 2 



(X s ,At/j)ds+ / (X s , 
Jo 



is a martingale with the same quadratic variation as for super-Brownian motion with variance pa- 
rameter a 2 . The Dawson-Watanabe process with killing rate (which we sometimes refer to as 
the standard Dawson-Watanabe process) is super-Brownian motion. Existence and distributiona l 
uniqueness of Dawson-Watana be processes in gen eral is asserted in IDawson and Perkins (1999) 
and proved, in various cases, in lDawsonl (1197 8l) and Evans and Perkins! (11995b . It is also proved in 
these articles that the law of a Dawson- Watanabe process with killing on a finite time interval is ab- 
solutely continuous with respect to that of a standard Dawson-Watanabe proce ss with the same vari 



ance p arameter, and that the likelihood ratio (Radon-Nikodym derivative) is (IDawson and Perkins 

(EHi)) 



(12) 



expj-y 9(t,x)dM(t,x)-^ J(X t ,9(t,-) 2 )dty 
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where dM(t , x) is t he ort hogonal martingale measure attached to the standard Dawson-Watanabe 
process (see Walshl (1986)). Absolute continuity implies that sample path properties are inherited: 
In particular, when d = 2, 3, if X t is a Dawson-Watanabe process with killing, then almost surely 
its occupation time process L t is absolutely continuous, with local time density L t (x) jointly con- 
tinuous in x and t. 

It is shown in lLallevid2007h that for the SIR- 1 epidemic in Z with village size N, the particle den- 
sity processes, suitably rescaled, converge as N — > oo to the density process of a standard Dawson- 
Watanabe process or a Dawson-Watanabe process with location-dependent killing, depending on 
whether the total number of initial infections is below a critical threshold or not. In dimensions 
d > 2, one cannot expect such a result to hold, because the Dawson-Watanabe process is almost 
surely singular with respect to the Lebesgue measure and therefore has no associated density pro- 
cess. However, as measure-valued processes, the SIR-d (d = 2, 3) epidemics, under suitable scaling, 
do converge, as the next theorem asserts. For the SIR-d model with village size 7Y, define 

(13) Xf {x) : = # infected particles at x at time i; 

(x) : = # recovered particles at x at time n = Xf {x). 



Theorem 2. Assume that d = 2 or 3, and suppose that for some a < 1/(3 — d/2) the initial 
configurations fj, N := Xq are such that 

(14) !Fj\[afi N /i with compact support, and 

(15) (fa N G Nat )(Vmx))/N a V- d M [(q aH * p)/a 2 ](x) e C(R 1+d ) 
where the second convergence is in D([0,oo);C b (R d )). Then 

(16) (j r N<*X N )NO: t ==> Xf 

where the limit process Xt is a Dawson-Watanabe process with initial configuration Xq = p, 
variance parameter a 2 , and killing rate 9. The killing rate depends on the value of a as follows: 

(i) if a < 1/(3 - d/2), then 6 = 0; and 

(ii) if a = 1/(3 - d/2), then 6 = L t (x), 

where Lt(x) is the local time density of the process Xt. The convergence => in (1161 ) is weak conver- 
gence relative to the Skorokhod topology on D([0, oo); M^(M d )). 

Theorem [2] will be proved in §3. 

Remark 5. Theorem [2] asserts that there is a critical threshold for the SIR— d epidemic in dimensions 
d = 2, 3: Below the threshold (when the sizes of the initially infected populations are N a ", 
where a* = 1/(3 — d/2) is the critical exponent) the effect of finite population size is not felt, and 
the epidemic looks much like its branching envelope. At the critical threshold, the finite-population 
effects begin to show, and the epidemic now looks like a branching random walk with location- 
dependent killing. 

Remark 6. The critical behavior of the SIR-d epidemics in dimensions d > 4 is considerably simpler. 
It can be shown that in all dimensions d > 5, the critical threshold is a = 1; in the four dimensional 
case, the critical threshold is still a = 1 but there will be logarithmic corrections - we thank Ed 
Perkins for pointing this out to us. 
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1.6. Notational conventions. Since the proof of Theorem |2] is based on likelihood ratio calcula- 
tions, we shall, at the risk of minor confusion, use the same letters X and R, with subscripts and/or 
superscripts, to denote particle counts and occupation counts for both branching random walks and 
the SIR-d epidemic processes (see equations (fT2l and (fT3l ) and for their continuous limits. Through- 
out the paper, we use the notation / x g to mean that the ratio f / g remains bounded away from 
and oo. Also, C, Ci, etc. denote generic constants whose values may change from line to line. The 
notation 5 x (y) is reserved for the Kronecker delta function. The notation Y n = op(f(n)) means 
that Y n /f(n) — > in probability; and Y n = Op(f(n)) means that the sequence \Y n \/f(n) is tight. 
Finally, we use a "local scoping rule" for notation: Any notation introduced in a proof is local to the 
proof, unless otherwise indicated. 

2. Local time for branching random walk m d = 2,3 

2.1. Estimates on transition probabilities. Recall that P n = (P n (x — y)) is the n— step transition 
probability kernel for the simple random walk on (with holding parameter l/(2d + 1)). For 
critical branching random walk, P n {x,y) is the expected number of particles at site y at time n 
given that the branching random walk is initiated by a single particle at site x. For this reason, sharp 
estimates on these transition probabilities will be of crucial importance in the proof of Theorem Q] 
We collect several useful estimates here. As the proofs are somewhat technical, we relegate them to 
the Appendix (section |4]below). Write 

$ n (x, y) = 4> n (x) + <j) n {y) where 

is the Gauss kernel in M. d . The first two results relate transition probabilities to the Gauss kernel. 

Lemma 2. For all sufficiently small (3 > there exists constant C = C(f3) > such that for all 
integers m, n > 1 and all x £ Z d , 

(17) P n (x) < C(p n {(3x) and 

(18) (P m * <M(J3x) < C<p m+n (Px/2). 

Furthermore, for each A > and each T > there exists C = C(A, T) > such that for all k 
sufficiently large and all \x\ < A\f~k, 

(19) ]T M(3x) < C P n(*)- 

n<kT n<kT 

Lemma 3. For all sufficiently small (3 > there exists constant C = C(/3) > such that for all 
integers n > 1 and all x,y € Z d , 

(20) \P n (x) - P n (y)\ < C ( liZjd A • Z> n (J3x,Py). 



n 



In particular, for all 7 < 1, 

(21) \P n {x) - P n (y)\ < C ( ) ■ $ n {l3x,l3y). 



, 7 

y ' 



Our arguments will also require the following estimates on the discretized Green kernel. 
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Lemma 4. For each 7 6 (0,2 — d/2), (3 > 0, h > 1, n G N, a«d s, y G de/wre 

(22) F n , h (x, y; /3) = F njh . 7 (x, y; /?) = £ ^ 7^ + p), /% + p)), 

|p|</i i<n 

where the first summation is over p G 7L d with \p\ < h. Then there exists C = C{pf, (3, d) < 00 such 
that for all n G N, hi, > 1, a«J all x,y G Z d , the following inequalities hold: 

(23) F n)fel (x,y;/3)-F n ^ 2 (x,y;^) < Cn 2 -^/ 2 F n ^i+^-i(^, 1/5 
^^-Pi(z) • [-^n-i^i^ - 2,y- z;/3) ■ F n _ iM (x -z,y- z;f3)] 

i<n z 

(24) < Cn 2 ^ d+2 ^' 2 ]T Y, ^ x + P)/ 2 > P(V + <°V 2 ) 

|p|</ii+/i 2 -l i<n 
< Cn 2-(d + 7)/2 Fn)hi+h2 _ i(X)y;/3/2) . 

(TVofe in ?/ze last term the (3 parameter is changed to (3/2.) 
Lemma 5. For each (3 > 0, h > 1, m, n G N, and x G define 

(25) Jm,n,h{x\l3) = Yl MP(x + p)), 

\p\<h m<l<m+n 

where the first summation is over p G Z d wz'f/i |/?| < /i. 77ie« f/iere exitfto C = C([3) > such that 
for all m,n G N, h\,h,2 > 1, awe? a/Z x G the following inequalities holds: 

(26) J m>n>hl (x; (3) ■ J m ,n,li2 
(27) 

• [ J m ,n-iM (x - z;P)- J mjn -i th2 (x - z;P)\ < Cn 2 ~ d/2 J m , nyhl+h2 -i(x; (3/2). 

i<n z 

2.2. Proof of Theorem [TJ For notational ease, we omit the superscript k in the arguments below: 
thus, we write X n (x) instead of X*(x), and R n (x) instead of R*(x). To prove the theorem it 
suffices to prove that (1) the sequence of random processes (Rkt(V~kx)/k 2 ~ d / 2 ) is tight in the space 
D([0, 00); Cb(M d )); and (2) that the only possible weak limit is the local time density process L t (x). 
The second of these is easy, given LemmaQ] This implies that for any test function -0 G C c (M. d ), 

-j-2^2 R kt(^kx)il;(x) 

X 

i<kt x 

f X s {il>)ds, 
Jo 
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where (X t ) is the super-Brownian motion started in configuration X = fi, run at speed a. Hence, 
any weak limit of the sequence (R kt (V~kx) / k 2 ~ d / 2 ) must be a density of the occupation measure 
for super-Brownian motion. On the other hand, by Remark[3]and Sugitani's Theorem, 

X s (ip)ds = / Lt{x)ip{x) dx. 



o 



It follows that L t (x) is the only possible weak limit. 

Thus, to prove Theorem Q] it suffices to prove that the sequence (R kt {\fkx) / k 2 ~ d l 2 ) is tight in 
the space D([0, oo); C{,(R. d )). In view of hypothesis ©, it is enough to prove the tightness of the 
re-centered sequence 

(28) Y k (t, x) := (R kt (Vkx) - ^ k G kt )(Vkx)) / k 2 ~ d / 2 . 



This we will accomplish by verifying a form of the Kolmogorov-Centsov criterion. According to 
this criterion, to prove tightness it suffices to prove that for each compact subset K of [0, oo) x M. d 
there exist constants C < oo, a > 0, and 5 > d + 1 such that for all pairs (s, a), (t, b) G K, 

(29) E\Y k (t,a) - Y k (t,b)\ a < C\a - b\ 5 and 

(30) E\Y k (t,a) - Y k (s,a)\ a < C\t - s\ 5 . 

The trick is to not work with moments directly, but instead, following the strategy of lSugitani dl989l) . 
to work with cumulants: 

Lemma 6 (Lemma 3.1 in ISugitanil d 19891) ). Let X be a random variable with moment generating 
function E exp(9X) = exp(^^ =1 n a n ). If for some integer N there exists r,b > such that 

K| < br n , forn < 2N, 

then there exists C = C(b, N) > such that 

EX 2N < Cr 2N . 



A. Cumulants. In the following discussion we use the notation (u, f) to denote the inner product 
of a function / and a measure v on Z d , and we let R n be the occupation measure of the nearest 
neighbor branching random walk with Poisson(l) offspring distribution. By the additivity and spa- 
tial homogeneity of the branching random walk, for any ip G C c (JL d ) and for each n > 1 there 
exists a function v n = v t G C c (Z, d ) such that for any (nonrandom) initial configuration 

E^exp((R n , V>)) = exp((/x, v n )). 

Note that v\ = ip. The assignment tp i— > v% is monotone in tp, but not in general linear. Setting 
H = 5 X and conditioning on the first generation, we obtain 

exp(z^ n+ i(x)) = Y^jQi ^ 2d + l X] exp ^( x + e ^ + Vn(x + e))j 

where {Qj} is the offspring distribution (in the case of interest, the Poisson distribution with mean 1) 
and the inner sum is over the 2d + 1 nearest neighbors e of the origin in Z d (recall that the origin is 
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included in this collection, since particles of the branching random walk can stay at the same sites 
as their parents). Observe that if the offspring distribution is Poisson(l), then 

(31) fn+i(s) = 2 <i + i ^2 eyL P(^( x + e ) + v n {x + e)) ~ 1- 

e 

Define the cumulants K^ n (x) = K^ n (x) in the usual way: 
By the arguments of the preceding paragraph, = tp and k^i = for all h > 2, and by (1311 ), 

X O k K h ,n+l(x) = eX P { X ° hK hA X + e ) + ^ + e) > - 1 ] • 

h>i e y [fc>i J y 

Consequently, 

(32) Kfc.n+^a;) = 2rf 1 X X ^7 X n^^."^ + e ) + ^(^i) " + e )i 

e m=l ' V m (h) i=l 

where V m (h) denotes the set of m— tuples (hi, h,2, ■ ■ ■ , /i m ) of positive integers whose sum is h, 
and 6\(-) is the Kronecker delta function. When h > 2, the m = 1 summand in (|32T ) equals 

l/(2d + 1) • ^^(i + e) = (Pi * Kh,n)(x), hence, 

Kh,n+i(x) = (Pi * K h>n )(x) + E n+ i(x), 

where 



j h ^ rn 

(33) 5 n+ i(g) = S n+ i(a;; ^) := M , ^ E^j X II^.™^ + e ) + M^i)^ + e )}- 

e m=2 m ' p m (ft) i=l 

Since k^i = for all h > 2, by iteration we then get 

n-l 

(34) k m (x) = J^(P * H n _;)(x). 

Z=0 

Consider now the special case V> = ^a,& := $a — &b where a, 6 G 7L d and ^ is the Kronecker 
delta function. Fix 0<7<2 — d/2 small, and let 

17 = 77(7) = 2- (d + 7)/2 > 0. 

Recall that in (1221 in Lemma|4]we defined F n ^(x, y; (3) for {3 > 0, /t > 1, n G N and x,y £ Z d as 

^(a, y; /3) = ^n,h ;7 (^, y; /?) = XX E72 + ^ ^ y + 

\p\<h l<n 
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Claim . For each h > 1 there exists < oo such that for all n G N and all x G 

(35) \K h , n (x)\ < C h \a - b\ h ^n^ h ^F rhh (a -x,b-x; 2'^ (5); 
moreover, for all h > 2, all n G N an d all x G 

(36) < Cfclo-fcl^n^- 1 )-^ 2 ]T ^ ^( 2 -(^ 1 )/3(a-x + /0 ), 2^^ 1 )/3(&-x+p)). 

|p|</i i<" 

In fact, when /i = 1, 

= |^(^,^)| = \G n (a - x) - G n {b - x)\ 
< C |a - 6|T ^ -^<&jG0(a " - *)) 

= C|a — 6| 7 F ni i(a — x,b — x; (3), 

where in the middle inequality we used assertion (|2TI) of Lemma|3] Furthermore, since t/>(a;) / if 

and only if x = a or 6, in which case |^(a;)| = 1 and 

inf n C\a — 6| 7 F ni i(a — x, b — x; (3) > 0, we get that for all n and all x, 

(37) |«i,»(aO + ^(x)| < C|a - 6| 7 F n ,i(a - x, b - x;0). 

Now suppose that the claim holds for 1, . . . , h — 1, and we want to prove the claim for /i. First note 
that in the definition (1331 of S n (x), only k^. for hi < h are involved, hence by induction, (1371 ) and 
relation (1231 we get that for all n and x, 

~ n (x)<C\a-b\^n^Y.i2 E 

e m=2V m (h) 

■ F nM (a-x + e,b-x + e; 2~^-^(3) ■ F n>h _ hl {a - x + e, b - x + e; 2- {h ~^0). 
The claims then follows from (l34l) and (l24l). 



B. Proof of d29l ). Suppose the initial configurations ^ satisfy the hypotheses of Theorem Q] For 
any a, 6 G M d , we want to estimate \ Yk(t, a) — Yj.(t, b)\ m . By Lemma|6l this can be done by 
setting ip = - and estimating \ (fi k , n h ,kt)\- B y OSJl. for all h > 2, 

|(/A«Mt>| < C A |>/fc(o- 6)|^ • (fct)^- 1 )- 7 / 2 - 

• ^ <//, ^ ^(2-(^ 1 )/?(v / fca + p - •), 2-^(3(Vkb + p- •)))• 

\p\<h Kkt 

By assumption (0]), we can find an A > such that supp(^ fc ) C B(0, Ay/k) for all k, where 5(0, r) 
represents the ball of radius r around 0. We therefore have that for all k sufficiently large, 

max<//, M2- (h - 1] P(x ~ ■))> = max £ £ / (2^ ( ^ 1) /?(x - y)) 
(38) <C^ W (y)^PKx-y) 

2/ 

< C^ 2 , 



BRW LOCAL TIME 13 



where in the first inequality we used ( fT9l , and the last inequality is due to the relative compactness 
of (/j, k Gkt)(Vk-)/k 2 ~ d / 2 assumed in ©. Hence when h > 2, 

\{» k ,Kh,kt)\ < Ck h ^ 2+r i( h - l ^l 2+2 - d ' 2 ■ f'^- 1 ^' 2 ■ \a - b\ h \ 
Plugging in 77 = 2 — (d + j)/2 gives us 

\(»\Kh,kt)\ < Ck<^ d ™ h • t (2-(«^))(h-l)-7/2 . | fl _ 6 |*7. 

Noting that {Y k (t, a) - Y k (t, b)) = 0, by LemmaH we get 

E » k \Y k (t, a) - Y k (t, b)\ 2h < Ct(2-(<i+7)/2)(2ft-l)-7/2 . _ 6 |2ft 7 _ 
By choosing h large such that 2/17 > d + 1 we obtain d29l ). 



C. Proof of (1301) . By the additivity and spatial homogeneity of the branching random walk, for any 
if) G C c (Z d ) and for all m, n G N there exists a function y n = u n ^ C c {'L d ) such that for any 
(nonrandom) initial configuration 

£ M exp((i? n+m - R m ,4>)) = exp((//,z/( njm ))). 

Letting [i = 5 X and conditioning on the first generation, we get 

exp(z/(„ im+ i)(a;)) = ^ exp((i?„ + i +m - i? m+ i,0)) 

where Q = {Qj}j>o denotes the offspring distribution. In case where the offspring distribution is 
Poisson(l), the equation above implies 



(39) 



"(n,m+l) ( X ) = 2 (i+ 1 ^ 6XP ^ {n ' m) ( X + e )) ~ L 



Define the cumulants K.h,(n,m) Dv 

^exp(e(i2^-i2 m ,V'))=exp((M > X) e V(n,m)>)' Vflei 



/i>i 

Then by <|35), 



e 



1. 



Therefore 



(40) Kf,,(n,m+i)(aO = ^TtZX^ Il K ^,Km)(^ + e), 

e i=l - Vi (h)j=l 

where T-^/i) denotes the set of z— tuples {hi, h%, ■ ■ ■ , hi) of positive integers whose sum is /i. The 
m = 1 summand in (|40l equals l/(2cZ + 1) • ^ e Kh( n ^{x + e) = (Pi * Khf n ,m))( x )> hence when 
/i > 2, 

Kh,(n,m.+l)(^) = (-Pi * Kh,(n,m)){x) + Hn,m+l(z), 
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where 

1 - 1 - 

(41) ~ n . m+ i(x) =E n>m+1 (x;h) := — ^ ^ - ^ JJ K ft . ){n!m) + e) 

e i=2 " v l {h)j=i 

By iteration we then get that for all h > 2, 

(42) K/i,(n,m)O0 = ^(Pj * H„ im _i)(x). 

For if) = l a for a G Z d , by C7]) 

«i,(n,*oO»0 = E s *(R n+m -R m ,ip) = Pi{a-x)<C <h(P(a-x)). 

m<l<m+n m<l<m+n 

Furthermore, similarly as in proving the claim in Part A, using Lemma|5]and (1421 we get that for all 
h, there exists > such that 

K, (n , m) (x)\ < C h n^^ £ * (2- (ft - 1} /?(« -* + p))- 

\p\<h m<l<m+n 



We are ready to verify (I30i Setting tp = and using (1381 ). we get 

_ Qj f; h(2-d/2) . t (2-d/2)(h-l) 

By Lemma [6] we then get 

E » h \Y k (t + s,a)- Y k (s, a)\ 2h < C ■ ^- W^-i) . 
So by choosing h large such that (2 — d/2)(2h — 1) > 1 + d we obtain (l30l ). 



□ 



Remark 7. For general offspring distributions Q = {Qj}, let f(x) = log(^ ■ QjX^) where x > 0. 
If the offspring distribution Q has an exponentially decaying tail, then f(x) can be expanded around 

x = 1 as /(x) = YT=i f [i \ l ){ x - i)V^ ! - Thus <GD turns into 

z, n+1 (x) = fj /W(l) ^ exp W(x + e) + i/ n (a; + e)) -l^j £\, 

and 

pVn + iW = f ^E ex P |% + e) + ^V(x + e)j -lj /V 

This enables us to express K/, in +i(x) in terms of ip(x + e) and Kh,n{ x + e ) similarly as in (l32l and 
in (l34l (note /^(l) = 1 because Q has mean 1), and prove the Kolmogorov-Centsov criterion 
for the spatial variable. Similarly one can verify the Kolmogorov-Centsov criterion for the time 
variable. 
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2.3. Sufficient conditions for Assumption ®. Now we state some conditions that imply (© and 
are easier to check. 

Proposition 1. Let d = 2 or 3. Suppose that the initial configurations fi k are such that J~kH k => /A 
a«<f satisfy 

(43) limsupmaxf/GfctHx)/*; 2 -^ 2 = 0. 

Then ([9]) holds. In particular, if any of the following assumptions is satisfied, then ((9]) holds. 

(i) In dimension 2, the maximal number of particles on a single site is bounded in k, i.e., 
sup fc max y p> k (y) < oo. 

( ii) In dimension 3, there exist C\,Ci > such that 

(44) C7 2 :=supmax fi k (y) < oo, 

where B(x, r) denotes the ball of radius r around x for any x and r > 0; that is, the number of 
particles in any ball of radius 2>Cik 1 /® is bounded in k. 

(Hi) In dimension 2, is such that fi k (y) is a decreasing function in \y\, and there exists a G (0, 2) 
such that 

(45) /(y) <C(V fc/(|y|2 + l)y\ Vy,k. 



Remark 8. This proposition is a natural analogue of Proposition 1 in lSugitanil (I1989I) . 

To prove Proposition [T] we will need the following result. 
Lemma 7. For any function ifi G C c (R. d ) and each integer k > 1, define 

^ k (x)= ^2ip(y/Vk)G kt (Vkx-y)/k, for x G Z d /Vk and t G Z/Jfe, 



and extend by linear interpolation elsewhere. Then 

(46) lim = [(q aH *i>)/o- 2 ](x), 

fc— »oo 

and the convergence is locally uniform in t and x. 

Proof. Pointwise convergence d46l ) follows from the local central limit theorem. To prove that 
the convergence is locally uniform, it suffices to show that the sequence of functions (^(x)) is 
relatively compact in C(R 1+d ). For this, we use the Ascoli-Arzela criterion. First, we show that the 
functions V?* (x) are uniformly bounded on any compact set in M} +d . Denote by M the maximum 
of Then 

< Yl Mv/Vk)\-Gkt(Vkx-y)/k 

( 47 ) <M Y J Gkt(Vkx-y)/k 

< Mt. 
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Next, we show they are equi-continuous. Fix e > 0, and set 5 = e/M. By ( |47T ), |*^| < e for all 
t < S; thus, 

<2e, Vx,y GM d ands,t < <S. 
On the other hand, by (|20l) . for all i > <5 and M 2 , 

k8<n<kt * z 
kS<n<k 

<2e + C5^ d ~^ 2 -\x-y\ 



(In the second inequality we used the fact that ^ < CVk d ; this holds because t/> is 

bounded and has compact support.) Finally, for all x and all <5 < s < t, 

<M ^ ^P n (V^x-^)/A; 

ks<n<kt z 
< M(t- s). 

□ 

Proof of Proposition^ For any ?/> £ C c (M d ), by Lemma|7J ^ (a;) converge to [(q , -2 t * ip)/a 2 \{x) 
in the local uniform topology. Therefore, 



k 
y 



(At, [(^ * if))/* 2 }) (-.- M fe (v^-)A =>■ M G M F (M rf )) 

= ([(9cr2t*A^)/^ 2 ]^}- 

On the other hand, if we can show that (fi k Gkt)(Vkx)/k 2 ~ d ^ 2 is relatively compact in C(]R 1+rf ), 
then for any limit F(t,x), 

Y,(v k G k t)(Vkx)/k 2 4>{x) = J2(n k G kt )(Vkx)/k 2 - d / 2 ■ ij){x) ■ 1/Vk d f F(t,x)^(x)dx. 

x x -> x 

Hence, {[{q a 2 t * fi)/a 2 ], if)} = J x F(t, x)ijj(x)dx, which implies that (1) the measure [(q a 2 t * p)/cf 2 ] 
has density, and (2) (fi k G kt ){Vkx)/k 2 -^ 2 converge to [(q a 2 t * fi)/a 2 ](x) in C(E 1+d ). 

Now we show that d43l ) implies that (p, k Gkt){Vkx) j 'k 2 ~ d l 2 is relatively compact in C(R 1+d ), 
by verifying the Ascoli-Arzela criterion. We first show that they are uniformly bounded on any 
compact set in M 1+ci . In fact, by (|43T ). there exists 5 > such that 

supm&x(fi k G k5 ){Vkx)/k 2 - d / 2 < 1; 

k x 
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moreover, for all i > 5 and all x, 

(fi k G kt )(Vkx)/k 2 ^ 2 

<i+ 2 S> fc (4P»(v'te-*)/fc 2 ~ <i/2 

kS<n<kt z 

<1 + ^ C7/n d / 2 • fc//fc 2 " d / 2 
fc<5<n<fct 

<C = C(t), 

where in the second inequality we used the fa cts that there exists C > such that for all n and all 
x £ Z d , P n (x) < C/n d / 2 rcf. ISpitzeJ dl976h . Proposition 6 on p72), and that the total number of 
particles Yl z ^ k ( z ) = 0(k). 

Next we show they are equi-continuous. In fact, for any e > 0, by (|43T ). there exists 5 > such 
that 

supm&x(n k G k s)(Vkx)/k 2 ~ d/2 < e; 

k x 

therefore, for all s, t < 5 and all x, y, 

snp\(fi k G ks )(Vkx)/k 2 - d / 2 - (fi k G kt )(Vky)/k 2 - d / 2 \ < 2e; 

k 

moreover, for all t > 5 and all x ^ y, by (|20l ), 

\{li k G ht ){Vkx)/k 2 - d l 2 - (fi k G kt )(Vky)/k 2 ~ d / 2 \ 



fc 2-d/2 

kS<n<kt z 



^ 2£+ CVk\x-y 



kS<n<kt z 



7^~, V™ 



fc<5<n<fct 



/n 

< 2e + C7r( d - 1 )/ 2 |x-y|; 
and for all x and all 5 < s < t, 

\(fi k G kt )(Vkx)/k 2 - d / 2 - (fi k G ks )(Vkx)/k 2 ^ 2 \ 

= ^i/2 E E^-») 

ks<n<kt z 
fcs<n<fci 



< 



Clog(i/s) < C(t-s)/S, i£d = 2; 

C(l/^-l/Vi)<C(t-s)/8 3/2 , if d = 3. 
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We have therefore proved that (l43T > implies the relative compactness of (fi k Gkt)(V~kx)/k 2 ~ d / 2 . 
Next we show that any of the conditions in (i)^(iii) implies (|43T >. 

(i) For all x £R 2 and all t > 0, 

(» k G kt )(x)/k = ^EE » k ( z ) p n{x ~ *) 

n<kt z 

<j E>=», 

n<kt 

therefore (|43T ) holds. 

(ii) In order to verify (l43l) . by (fT71) . it suffices to show that 

(48) limsupmax fi k ly)4> n ((3(x - y))/Vk = 0. 

t— >0 ju a: 

Claim . There exists C3 > swc/i that for all k, n and all pairs x, y E "L d with \x — y\ > C\k 1 ^, 
Vk~4>n(P(x - v)) < c 3 4>n{0{x- Z )). 

In fact, the above inequality is equivalent to 

Vk<C 3 £ e,J ^ X -^-^- Z ^ ),y\x-y\> Cl k^. 

But this holds trivially since when x ^ B(y,C\k 1 ^ 6 ), there is a positive proportion of integer 
points z in the ball B(y, Cifc 1 / 6 ) such that |x — y\ > \x — z\, and the proportion does not depend on 
k, x, y. Now let us estimate ^2 n<kt ^2 y n k (y)4> n (P(x — y))/ \fk: For any fixed k and x, this sum 
can be written as the sum of the following two terms: 

/:= E E M fc (i/)&»(/?(s-t/))/Vfc, 

and 

//:= E E Ay)KiP(x-y))/Vk. 

n<kt [y-x^dk 1 / 6 

As to term I, we have 

1 ^ E E Ai/) • C/n 3 / 2 /^ 

< C ? ■ C/n^/Vk 

n<kt 

< C/Vk, 
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where in the second inequality we used (l44l ). And by the claim and (l44l ). 

11 ^E E Av)/^-C 3 £ <t> n {P(x-z))/yfk 

n<kt \y- x \>C\k 1 / 6 zeB{y,Cik 1 l 6 ) 
n<kt 2 yeBiZtdk 1 / 6 ) 

< E c a 

n</c* 

< Ct. 



Therefore (08]) holds. 

(iii) In order to verify (|43T ). by (fT71) . it suffices to show that 

(49) limsupmax fj, h (y)4> n (/3(x - y))/k = 0. 

t— >0 u x — ' * — ' 

n<kt y 

By assumption, n k {y) is a decreasing function of |y|; so is 4> n {[3y). Therefore, by Lemma [8] below, 
the last term is bounded by ^2 n<kt iJL k (y)(j) n (Py)/k, which, by assumption (l45l) . can be further 
bounded by 



n<kt y \ V |y| / 

ht n ■ k Jxi 



n<kt 



n<fct 

< C f s~ a/2 ds 



x\~ a e 2 "/ fe dx 





= 0(t), 

where the third inequality and the last equation hold because a < 2 by assumption. □ 

Remark 9. In dimension 2, if the assumption in (iii) is satisfied, then the radius of the support of n k 
will be of order y/k. This is because we need J2 y ^ k {y) = 0(k), hence for some C > 0, 

E (Vk/(\y\ 2 + l)) a >Ck, 

»/gSupp(/i fe ) 

i.e., 

E (vVW+^) a > ck 1 -^ 2 . 

j/esupp(/i fe ) 
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But for any r, 



(i/vwn) a = o( [ (1/VW+1Y dy] 



( J 

2-a 



= 0(r 

so in order that 0(r 2 ' a ) > Ck 1 ~ a / 2 , we need r = Oik 1 / 2 ). 

Remark 10. In dimension 3, if fi k is such that /i k (y) is a decreasing function in \y\, and there exists 
a G (0, 2) such that 

(50) fi k (y) <C^k/(\y\ 2 + l)Y/Vk, Vy,k, 



then using the similar proof as in (iii) we can show that (1431 holds. But in fact, (1501 ) can never be 
satisfied. The reason is that in order fi k (y) > 1, we must have that for some C > 0, 

l<c(V£;/(|2/| 2 + l)) a /v^ 

which implies that 

\y\ < Ck 1 ^ 1 ^ = o(k 1 ' 4 ). 
On the other hand, we need fi k (y) = O(k), therefore for some C > 0, 

£ (Vk/(\y\ 2 + l)) a /Vk>Ck, 

|»|=o(fcV*) 

or 

(5D e (vvwn) a >ck^/ 2 . 

|j,|=o(fcl/4) 



However, 



1,1=^/4) v ; VA,i=o(^ 4) v ; / 



ro(fcV*) „ N 

O I / (1/Vr 2 + 1) -r 2 dr 



= (jfe(3-«)/4) = (jb(3-a)/2) 5 

contradiction with (T5TT) . 

Lemma 8. Suppose that f and g are two nonnegative functions on Z d , and f has compact support. 
Suppose further that both f(x) and g{x) are decreasing functions in \x\, then 

(52) ^f(y)g(x-y) <^f(y)g(y), Vx. 
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Proof. Since f(x) is a decreasing functions in \x\ and has compact support, we can enumerate its 
positive values, say a\ > . . . > a n > 0. We can also enumerate the values of g, say b\ > . . . b n > 
.... To show d52~l ), it then suffices to show that 

n n 

sup y^dkbjt, = yjqfc&fc. 
h,-,in k=1 k=l 

But this is easily seen to be true. □ 
3. Proof of TheoremEJ Spatial epidemics in Dimensions d = 2,3 



3.1. Strategy. The strategy is the same as that used by iLallevi (2007) in the 1— dimensional case 



Since the law of the SIR-d epidemic with village size N is absolutely continuous relative to that of 
its branching envelope, and since the branching envelopes converge weakly, after renormalization, 
to super-Brownian motion, it suffices to prove that the likelihood ratios converge weakly to the like- 
lihood ratio ([121 of the appropriate Dawson-Watanabe process relative to super-Brownian motion. 
The one- and higher-dimensional cases differ only in the behavior of the occupation statistics that 
enter into the likelihood ratios. 

3.2. Modified SIR-d epidemic. As in the one-dimensional case, it is technically easier to work 
with the likelihood ratio for a modification of the SIR- d epidemic. Recall that (a) when an infected 
individual attempts to infect a recovered individual in an SIR epidemic, the attempt fails; and (b) 
when two (or more) infected individuals simultaneously attempt to infect the same susceptible in- 
dividual, all but one of the attempts fail. Call an occurrence of type (a) an errant attempt, and an 
occurrence of type (b) a collision. In the modified SIR epidemic, collisions are not allowed, and 
there can be at most one errant attempt at any site/time. A formal specification of the modified SIR 
epidemic uses a variation of the standard coupling described in section [T72l as follows: 

Modified Standard Coupling: Particles are colored red or blue; red particles represent infected 
individuals in the modified SIR epidemic. Each particle produces a random number of offspring, 
according to the Poisson(l) distribution, which then randomly move to neighboring sites. Once 
situated, these offspring are assigned colors according to the following rules: 

(A) Offspring of blue particles are blue; offspring of red particles may be either red or blue. 

(B) At any site/time (x, t) there is at most one blue offspring of a red parent. 

(C) Given that at site x and time t there are y offspring of red parents, the conditional probability 

k n(v) = KN,t,x(y) that one of them is blue is 

(53) K N (y) = {yR/N} A 1, where 

(54) R = R?(x) = Y,Ys N (x) 

and Y t N (x) is the number of red particles at site x in generation t. (Thus, R = R^ {x) is the 
number of recovered individuals at site x at time t.) The red particle process is the modified 
SIR epidemic. 

Proposition 2. For each N > 1, versions of the SIR epidemic and the modified SIR epidemic 
can be constructed on a common probability space in such a way that ( i) the initial configurations 
p, N of infected individuals are identical, and satisfy the hypothesis (1141 ) of Theorem [2} and ( ii) the 
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discrepancy D t (x) between the two processes at site x and time t (that is, the absolute difference in 
number of infected individuals) satisfies 

(55) maxAW = o P (N a ). 

t,x 

This implies that after Feller-rescaling, the SIR-d epidemic and the modified SIR-d epidemic are 
indistinguishable. Consequently, to prove Theorem [2] it suffices to prove the corresponding result 
for the modified epidemic. 



Proof of Proposition^ This is essentially the same as the proof of Proposition 3 in lLalleyl (|2007l) . 
except that different estimates for the numbers of collisions and errant infection attempts in the 
SIR-d epidemic are necessary. These are given in Lemma [9] below. □ 

Lemma 9. For each pair (n, x) E N x Z d , let (x) and (x) be the number of collisions and 
the number of errant infection attempts, respectively, at site x and time n in the SIR-c? epidemic 
with village size N. Assume that the hypotheses (114l) - (115l) of Theorem [2] are satisfied, for some 
a < 1/(3- d/2). Then 

(56) £ Y, K (x) + « (*) - 1)+} = o P {N a ). 

n x 

The proof of this lemma makes use of the following result. 



Lemma 10. [Proposition 28 i mLallev and Zhens\ n2007\) I Denote by U n (x) the number of particles 
at x at time n of a BRW started by one particle at the origin, then 

n-l 

(57) EU n {xf = P n {x) +t 2 Y,Y. p i( z ) p n-i( x - *), 



i=0 

r 2 



where a is the variance of the offspring distribution. 

Proof of Lemma® Since the life length of the process is O p (N a ), it suffices to show that for any 

t > 0, 

n<N a t x 

Consider first the number (x) of collisions at site x and time n. For any susceptible individual 
rj, a collision occurs at i] if and only if there is some pair £, £ of infected individuals at neighbor- 
ing sites that simultaneously attempt to infect rj. Therefore given the evolution up to time n, the 
conditional expectation of T^ +1 (x) is bounded by C(J2 e ( x + e )) 2 /N- We want to show that 

(58) £ £pe(*)) 2 /iV = 0p (N a ). 

n<N a t x 

By the dominance of BRW over SIR epidemic, if we denote by U n (x) the number of particles at x 
at time n of a BRW started by one particle at the origin, and Xi, i = 1,2,..., the position of the 
initial particles of our epidemic model, then 

E{X^{x)f < Y E Un(x-Xi) 2 + 2Y / Pn(x-X i )P n (x-X j ), 
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which, in dimension 3, by Lemma [TOl can be bounded by C^ i P n (x — + 2^ P n (x — 
Xi)P n (x — Xj). Therefore 

£ Y, E ( X n(x)) 2 <C E (CiV a + C7iV 2 7^ 3 ) 

(59) n<N a t x n<N a t 

= 0{N 2a ), 

which is o(iV x N a ) since a < 2/3. In dimension 2, again by Lemma [TOl EU n (x) 2 < C(l + 
log n)P n (x), therefore 

^E(X^(x)) 2 <C (C(l + \ogn)N a + CN 2a /n) 

(60) n<N a t x n<N a t 

= 0(N 2a log AT), 

which is also o(iV x N a ) since a < 1/2. 

Now consider the number (x) of errant infection attempts at site x and time n. In order that 
there be more than one errant attempt, either (i) two or more infected individuals must simultane- 
ously try to infect a recovered individual, or (ii) infected individuals must attempt to infect more 
than one recovered individual. The number of occurences of type (i) during the course of the epi- 
demic is op(N a ), by the same argument that proved (T58T ). Thus, it suffices to bound the number 
of errant attempts of type (ii). This is bounded by the number B N (x) of pairs g, g' of recovered 
individuals at site x and time n that are subject to simultaneous infection attempts. Clearly, 

£,e C,q' 

where the sums are over all pairs ((£, g), (£, g')) in which g, g' are recovered individuals at site x 
and time n and £, £ are infected individuals at neighboring sites, and Z^ g and are independent 
Bernoulli(l/((2d + 1)N)). Hence, 

E(B» +1 (x) \G n )<C (j^X^x + e)) (R%(x)/N) 2 . 

By Theorem [J and the dominance of BRW over SIR epidemic, for all e > 0, there exists C > 
such that with probability > 1 — e, 

max< ai (i) < C N a ^ d ' 2 \ 

X 

Note further that 

E E^(E^(* + e )) < c E E^^)) 2 , 

n<N a t x \ e / n<N a t x 

which, by (l60l and (|59l , is bounded by CN 2a log iV in dimension 2 and CN 2a in dimension 3. 
Therefore, by enlarging C if necessary we have that with probability > 1 — 2e, the following holds: 

2 

N {x)/N) 2 < CN 2a ^ 2 ~ d ^N 2a log N/N 2 = o{N° ,. 

□ 



E E(E^(* + e))V 

n<N a t x \ e / 
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3.3. Convergence of likelihood ratios. In view of Proposition |2j to prove Theorem |2] it suffices 
to prove the corresponding result for the modified SIR epidemic defined in 93.21 For this, we shall 
analyze likelihood ratios. Denote by Q N the law of the modified SIR epidemic, and by P N the 
law of the branching envelope. Recall (cf. the modified standard coupling) that in the modified SIR 
process there can be at most one errant infection attempt, and no collisions, at any site/time x, t. 
Given the evolution of the process up to time t — 1, infection attempts at site x and time t are made 
according to the same law as are offspring in the branching envelope; the conditional probability 
that one of the attempts is errant is ftjv(y) (see equation (|53T>). Consequently, the likelihood ratio 
dQ N /dP N at the sample evolution X N := {X?(x)} x>t is 

«n d Q N IT TT P(y\^(l ~ K x(y)) +p(y + l|^)^Ar(y + 1) 



dP N 11 11 p(y\X) 



where 



y = x t N (x), 

A = Af(x) = J2x t N ^ 1 (x + e)/(2d+ 1), and 

e 

p(k | A) = \ k e~ x /k\ 
By the same calculation as in Lallevl d2007h . equation (53), this can be rewritten as 



(62) ^ = (1 + e N ) exp i - EE Af (x)fif (s) - i £ E *W } , 
where 

Af(x) := (Zf (*) - Af (*))/JV« 
e f{x):=R^(x)/N 1 - a ; and 
£tv = o p (l) under P^. 

That the error term is o p (l) follows by an argument nearly identical to the proof of Lemma [9] 

Observe that under P N , the increments (in t) of the first sum in the exponential constitute a 
martingale difference sequence. Furthermore, the quantities A^(x) in equation (|62~1) are the atoms 
of t he orthogona l martingale measures M N associated with the branching rando m walks X N 



See lLalleyl (120070 for the analogous representation in the one-dimensional case, and IWalshl (119860 
for background on stochastic integration against orthogonal martingale measures. The martingale 
measures M N can be defined by their actions on test functions tp G C^°(M. d ). Write (p, tp) for 
the integral of ip against a finite Borel measure p on W 1 , and for the Feller- Watanabe rescaling 
operator ©; then 

M t N (V) = (TncX^i/j) - (F Na xF,ip) - [\F Na X% aa ,A Na il>)d8 

Jo 

where A\~ is the difference operator 

A k il>(x) = E^ + e /^) " ( 2d + WW ) I i( 2d + 1 ) A: ~ 1 ] ' 
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The first sum in the exponential of equation d62l can be expressed as a stochastic integral against 
the orthogonal martingale measure M N : 

(63) ]T ]T A?(x)q?(x) = [ [ 6 N (t,x)M N (dt,dx), 

t>i xez d 

where 

qN u _\ _ r,N /. / ATn„\ I jvrl-a 



'(t,x) =R% at (VN°x)/N 1 

Proposition 3. Let X be the Dawson-Watanabe process with initial configuration fi and variance 
parameter a 2 , and let M(dt, dx) and Lt(x) be the associated orthogonal martingale measure and 
local time density process. Then under P , given the hypotheses of Theorem\2\ as N — > oo, 

(F N °X N , 9 N , M N ) (X, 0, M) if a < 1/(3 - d/2) and 
(F N <*X N , 6 n , M n ) =^ (X, L, M) if a = 1/(3 - d/2). 

Proof. Given the weak convergence of the second margin 9 N , th e joint convergence of the triple 



follows by the same argument as in Proposition 4 of lLalleyl (|2007h . The asymptotic behavior of the 
processes 9 N follows from Theorem Q] □ 

Corollary 1. If a < 1/(3 - d/2) then under P N , as N oo, 

dQ N 

(64) tt — ► 1 in probability 

dP" 

provided that the hypotheses of Theorem\2\on the initial configurations are satisfied. 

Proof. Proposition [3] implies that the sums (l63l converge to zero in probability as N — > oo. That 
the second sum in the likelihood ratio (l62l) also converges to zero in probab ility follows b y the same 



argument as in the one-dimensional case (see the proof of equation (60) in lLalleyl ( 2007 )). □ 



Proof of Theorem^ Corollary Q] implies that the modified SIR epidemics have the same scaling 
limit as their branching envelopes when a < 1/(3— d/2). Thus, to complete the proof of Theorem[2l 
it suffices to prove the assertion (fT6l) when a = 1/(3 — d/2). For this, it suffices to show that the two 
sums in the exponential of equation (l62l converge to the corresponding integrals in the exponential 
of equation (fT2l . The convergence of the first sum follows fro m Proposition |3] and the representation 
(l63l . By the same argument as in the proof of Corollary 4 of lLallevid2007h . 



(65) J2^Z A n( x )en(x) = / e N (t,x)M N (dtdx) =► / L t (x)M(dtdx). 
The convergence of the second sum 

(66) A N := £ £ ^(xfe^ixf =► J (X t , {L t f) dt 



follows by an ai - gument similar to the proof of equation (60) in lLallevid2007h . The idea is that if one 



substitutes the conditional expectation (x)/N 2a = E(A%(x) 2 \ Q n -i) for the quantity (x) 2 
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in the sum (|66l ) then the modified sum converges; in particular, by Theorem Q] and Watanabe's 
theorem, 

n x 

= ^EEE^-i^+ e )/( 2 ' l + 1 )]/ r >< [^(*w a(2 - d/2) ] 2 

n x e 



(X t ,(L t ) 2 )dt, 

where the second equation holds because a = 1/(3 — d/2). Therefore, it suffices to show that 
replacing (x) 2 by its conditional expectation has an asymptotically negligible effect on the sum, 
that is, 

A N - B N = oJl). 



By a simple variance calculation (see lLalleyl (|2007h for the one-dimensional case), this reduces to 
proving that 

(67) E J2( X n(x)f/N 4a x [R»( X )/N^-^]* = 0,(1). 

n x 

In fact, for all e > 0, there exists C > such that with probability > 1 — e, 

max< t (x) < CN^ 2 - d ' 2 \ 

X 

Note further that 

E E^E^(*+ e )] 2 ^ E E^(*)) 2 ' 

n<N a t x e n<N a t x 

which, by (l60l ) and (|59l ), is bounded by CN 2a log N in dimension 2 and CN 2a in dimension 3. 
Therefore, by enlarging C if necessary we have that with probability > 1 — 2e, the following holds: 

E E( A « ( x )fl NAa x i R n {x)/N al ?- d M] A < CN 2a logN/N 4a = o(l). 

n<N a t x 

□ 

4. Appendix: Proofs of Lemmas (HE] 

4.1. Proofs of Lemmas |2]-[3j The strategy is to consider the regions |x| < (2Lnlogn) 1//2 and 
\x\ > (Lnlogn) 1 / 2 separately. We begin with the unbounded region. By Hoeffding's inequality, 
since the increments of S n are no larger than 1 in modulus, 

P n {x) < P(\S n \ > \x\) < 2dexp(-\x\ 2 /(2dn)). 

Now for < P < l/Vd and L = L(f3) sufficiently large, 

exp(— |x| /(2dn)) < exp(— /3 2 |x| 2 /(2n))/n( d+1 )/ 2 , V \x\ > \J Ln log n. 

Thus, 

(68) P n (x) < Cexp(-p 2 \x\ 2 /(2nj)/n id+1)/2 = C<f>n{px)/y/n, V|a?| > ^Lnlogn, 



BRW LOCAL TIME 



27 



and 

(69) \P n (x) - P n (y)\ < C$ n (px,Py)/^ 



<c(^-^-Al\ $ n (/3x,Py), V\x\,\v\ > y/Lnlogn. 



This proves inequalities (fTTT ) and (l20l) for x and y outside the ball of radius (Ln log n) 1 / 2 . 

To deal with the region \x\ < (2 Ln log re) 1 / 2 we shall use the following crude estimate, valid for 
all points x G 7L d (Theorem 2.3.5 in lLawler and Limic ( 2007 )): 

\P n (x) - a- d <j> n (x/a)\ < C/(V^ d ■ n). 

For p = f3(L) > sufficiently small, 

<M/to) > l/{Vn-n d/2 ), V |x| < y/2Lnlogn; 

consequently, 



(70) \P n (x) - a- d 4> n (x/a)\ < C</) n (Px)/y/n, V|x| < y/2Ln\ogn. 

This obviously implies ( fTTT ) for x in the region |x| < {2Ln\ogn) l l 2 , and hence, together with the 
argument of the preceding paragraph, completes the proof of ( fTTT ). 

Similar arguments can be used to establish inequality d20l for points x and y in the ball of radius 
(Ln log re) 1 / 2 centered at the origin. First, it is easily seen that for sufficiently small f3 > 0, 

-0n(tf)| <C((|x-y|/^)Al)<I> n (/3rE,/3y), V x, y G M d 

Hence, by d70l ), d20l holds for for x and y in the ball of radius (Lnlogre) 1 / 2 . Therefore, since 
the choice of L is arbitrary, to complete the proof of (l20l it suffices to consider the case where 
|x| < (Lnlogre) 1 / 2 and \y\ > (2Lnlogre) 1 / 2 . In this case, choose a point z in the annulus \z\ G 
((Lrelog n) 1//2 , (2Ln log re) 1 / 2 ) such that |x — z\ + \z — y\ < 2\x — y\. Using the fact that (l20l) holds 
for each of the pairs x, z and z, y, we have 

|P n (x) - P n (z)| < C ((|x - z\/y/n) A 1) $ n (J3x,Pz) < 2C ((|x - z\/y/n) A l) $ n (Px,0y) 

and 

|P n (» - P n (y)| < C (\z - y\/yfc A 1) $ n (Pz,py) < C ((\z - y\/yfc) A l) $ n (Px,Py). 
Consequently, 

\P n (x)-P n (y)\<C((\x-y\/y/n)Al)-$ n (J3x,Py). 
This completes the proof of d20l . Inequality (1211 ) obviously follows from d20l . 

Proof of (fl9l when d = 3. The following argument works for all d > 3. Firstly, ^n<fcT 4>n{Px) 
is bounded by X^^L i 4>n(P x)- This is a decreasing function in |x|; moreover, by Lemma 4.3.2 
in lLawler and Limid d2007h . it equals Ci/\x\ d ~ 2 + 0(l/|x| d+2 ) as |x| -> oo for some C\ > 0. 



Secondly, for all k sufficiently large and all |x| < Ay/k, 

n<fcT kT/2<n<kT kT/2<n<kT 



28 



STEVEN P. LALLEY AND XINGHUA ZHENG 



note further that 

n>kT n>kT 

therefore there exists 5 > such that all k sufficiently large and all \x\ < A\fk, 



n<kT n=l 



For the nearest neighbor random walk, Yl^Li p n(%) equals C2/\x\ d 2 + 0(l/\x\ d+2 ) as \x\ — > oo 
for some C2 > 0. Relation (fl9l ) follows. 

Proof of (fl9l) when d = 2. In this case, one can deduce from the proof of Theorem 4.4.3 in 



Lawler and Limid (120071) that there exist Q > 0, i = 1, 2, 3, 4 such that for all \x\ < A\fk, 

Y MM ~ C 1 +C 2 log(kT/\x\ 2 ), 

n<kT 

and 

^ PnOz) ~C 3 + C 4 log(kT/\x\ 2 ). 

n<kT 

02J follows. 

To complete the proof of Lemma|2l it remains to prove inequality (fT8T ). 

Proof of (HU). By dTTb . it suffices to show that there exists C > such that for all x £ Z d and all 

i,j G N, 

(71) X>i(/Wj(/?(* - y)) < C0 i+i (^/2). 

v 

For all y G Z d , Let be the cube centered at y with length 1, and define 



Uv)= / (27ri/p 2 )- d / 2 exp(-p 2 \z\ 2 /(2i))dz. 
JzeQ(y) 

Then there exists C > such that for all i and all x, 

< cfa{x). 



Therefore to show (1711 ). it suffices to show that there exists C > such that for all x G Z and all 

i,j G N, 

(72) ]T ^(y)^(x - y) < Ccj) l+j ((3x/2). 

y 

Note (4>i(-)) is the probability mass function of the random variable [Aj], where Aj ~ N(0, i/ft-Id), 
and for any 2 G M d \ U y dQ(y), [z] is the unique y such that z G Q(y) (Aj takes values on U y dQ(y) 
with probability 0, so [Aj] is well defined a.s.). Hence J2 y 4'i{y)4'j(' ~ u) i s the probability mass 
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function of [Aj] + [Aj] with A, and Aj being independent. Since [Aj] + [Aj] differs from Aj + Aj 
by at most 2, 

Y^Uy)~H*-v) < [ (27T(i+j)/(3 2 )- d ^e W (-(3 2 \z\ 2 /m + mdz. 

y J\Z-X\<2 

It is easy to see that the last term can be bounded by C4>i + j((3x/2) for some C independent of i, j 
and x. 

□ 

4.2. Proof of Lemma |4j By the monotonicity of 4> n (x) in \x\, for all integers m,l > 1 and all 
x, y G M. d we have 

(/>m(x)(/)i(y) < <j> m (x)<j)i(x) + (j) m (y)(j)i{y) 

< C{ml)~ d/i (<t) m i/ im+l) {x) + 4>ml/(m+l){y))- 

Now note that for any t > and any x, 

Mx) = (2nty d / 2 exp(-|x| 2 /(2i)) < 2 d ' 2 ■ fa{x), 

and when t > 1, 

Mx) < <P\t](x) ■ (\t]/t) d/2 < 2 d/2 <f> m (x), 

where \t] stands for the smallest integer bigger than or equal to t. Further note that when m, I > 1, 
ml/(m + I) > 1/2. Using the two inequalities above we then get 

(73) 4>m{x)My) < C{miy d,i {4>\2ml/(m+l)-\ 0) + <i>\2rrd/{m+tj\ (v)) ■ 

and 

$ m (x,y)®i(u,v) = <j) m (x)(j)i(u) + (j> m (x)(j)i(v) + (t>m(.y)(f>l( U ) + <t>m{y)(l>l{v) 
< C(miy d/4 ($(- 2 mi/(m+01 ( x > y) + $ r2mV(^+01 ( U ' v )) 

Therefore for all hi,h 2 > 1, 

^n,/n (2;, y; 0)F n ,h 2 (x, y; (3) 

= E E E ^(ml)^ /2 «f m (/?(x + Pi), + Pi)) • $1(^(3: + p 2 ), + ft)) 
|pi|<hi |P2|</i2 m<n Kn 

,74) * E E EE cm)-* 4 -' 2 

|pi|<hl |p2|</i2 m<n Kn 

■ {®\2ml/(m+l)-} + Pl),P(y + Pi)) + *r2ml/(m+i)l + Pn)>P(V + »))} 

^ E • E E c ( m/ )"' i/4 " 7/2$ r2 m //(m. + oi(/ 3 ^ + - )'/ 3 (y + / ))- 

|p|<hi+/i2— 1 m<nl<n 
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Observe that when m, I G [1, n), ml/(m + Z) G [1/2, n/2), hence the last term is bounded by 

|p|<hi+h 2 -l Kn 

< Cn 2-(d +7) /2 . ^2 r^ 2 $Mx + p),Kv + p)) 

\p\<hi+h2-l Kn 
=Cn 2-(d + 7)/2 Fn ^ i+/i2 _ i(X)y;/?); 

i.e., (g3) holds. 

We now prove <[24]>. By (|74b . 

^ ^2 p i ( z ) • [-Fn-t./n (x- z,y - z; f3)F n _ iyh2 (x - z,y - z; (3)) 

i<n z 

sEE^w E E E 

i<n z \p\<hi+h,2—l m<n-i Kn—i 

■ C(m/r d / 4 -^ 2 $ r2mV(m+/)1 (J3(x -z + p),P(y-z + p)) 

<EE«^ H/a E • E 

m<n i<n |p|</ii+/i2— 1 i<min(n— m,n— 1) 

■ Yl P i(. Z )®\2ml/(m+l)] (/3(x - Z + p), /3(y - Z + p)). 

Using relation (fT8T ) and noting that \2ml / (m + 1)1 < max(m, Z), we can further bound the last term 
by 

Cn 2-(d + 2 7 )/2 £ . £ M(3{x + p)j% 0{y + p)/2) 

|p|</ll+/l2 — 1 *<W 

<Cn^ d *iy 2 F nM+ha -. x {xMP/2). 



□ 



4.3. Proof of LemmaHJ For all hi, hi > 1, all ieZ d and all integers m, n > 1, by d73l . 

Jm,n,hi {x, P)Jm,n,h2 (.X, /?) 

= E E ■ E 0/ 1 (/?(x+pi))0/ 2 (/3(x + P2 )) 



|pi|</ii \p2\<h.2 m<h,h<m+n 

(75) 



f E E ■ E 



|pl|</il \p2\<h-2 m<h,l2<rn+n 

■ (hh)~ d/A {0raii a /(h+fa)l 0% + Pi)) + ^r«ifa/(h+ia)l + P*))} 

ih/(h+h)] (P( x + p)). 

\p\<hi+h2— 1 m<h,l2<m+n 

Note that when Zi, Z 2 G [m, m + n), |~2ZiZ 2 /(Zi + Z 2 )] G [m, m + n), hence the last term is bounded 
by 

Cn 2-d/2 J- ■ MP(x + p)) = Cn 2 - d / 2 J m , nM+h2 _ 1 (x;f3), 

\p\<h\+h2 — 1 m<l<m+n 
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i.e., ((261) holds. 

We now prove (|27l) . By (1751) . 

^ ^ ^ ^ Pi( z ) ' {Jm,n—i,hi P)Jm,n—i,h2 { x 2 > P)) 

i<n z 

<EE^) E ■ E r2/i/ 2 /(ii+«2)i - z + p)) 

i<n z |p|<ft.i+ft.2 — 1 m<h,l2<m+n—i 

m<li,l2<m+n \p\<hi+h,2 — 1 i<min(m+n— /l ,m.+n— 12) 

Using relation (fTST ) and noting that [2^/2/(^1 + ^2)] G [min(/i, /2), max(/i, /2)], we can further 
bound the last term by 

Cn 2-d/2 £ . £ fa{p(x + P )/2) = Cn 2 - d ^J m , nM+h2 _ 1 (x;p/2). 

\p\<hi+h2-l m<i<m+n 

□ 
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